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ABSTRACT 

We consider the stability of an accretion disk wind to cloud formation when subject to a central radiation force. 
For a vertical launch velocity profile that is Keplerian or flatter and the presence of a significant radiation pressure, 
the wind flow streamlines cross in a conical layer. We argue that such regions are highly unstable, and are natural 
sites for supersonic turbulence and, consequently, density compressions. We suggest that combined with thermal 
instability these will all conspire to produce clouds. Such clouds can exist in dynamical equilibrium, constantly 
dissipating and reforming. As long as there is an inner truncation radius to the wind, our model emerges with a 
biconical structure similar to that inferred by Elvis (2000) for the broad line region (BLR) of active galactic nuclei 
(AGN). Our results may also apply to other disk-wind systems. 

Subject headings: hydrodynamics — instabilities — turbulence — galaxies:active — (galaxies:) quasars: 
emission lines — galaxies: Seyfert — (galaxies:) cooling flows 



1. INTRODUCTION 

Non-spherical outflows are a ubiquitous feature of AGN. In 
this Letter we discuss the stability of such an outflow launched 
normally from the disk with a decaying power law velocity pro- 
file. The motivation for this work is to understand the condi- 
tions under which the flow may be unstable to cloud formation, 
and to infer the structure of the region in which such clouds 
may reside. Elvis (2000) has argued that data require AGN 
BLR clouds to reside in a narrow biconical structure. Our work 
herein supports this possibilty. 

A variety of launching mechanisms for AGN accretion disk 
outflows have been studied: radiatively accelerated outflows 
(Arav et al. 1994), hydrodynamic line-driven winds (Proga et 
al. 1999, 2000), hydromagnetic disk winds (Konigl & Kartje 
1994), (Pelletier & Pudritz 1992), thermal wind-type out- 
flows in low luminosity AGNs (Pietrini & Torricelli-Ciamponi 
2000), and others. In this work we parameterize the disk wind 
independent of the launching mechanism. 

A disk wind outflow can become unstable from linear and 
nonlinear perturbations in velocity, density, temperature, ion- 
ization state, etc. These can lead to cloud formation. The cloud 
model is a leading, though not unanimously accepted, paradigm 
for the structure of the AGN BLR regions: models exist with 
(Elvis 2000; Urry & Padovani 1995) and without (Murray & 
Chiang 1995; Murray et al. 1995) clouds. The main reason for 
the lack of agreement is the uncertainty about cloud survival 
in AGN enviroments (Mathews & Ferland 1987). They may 
be unstable to evaporation (e.g., see (Pier & Voit 1995)), or 
be insufficiently supported by pressure and be dynamically un- 
stable to shredding (Mathews & Ferland 1987; Poludnenko et 
al. 2002). Among the suggested solutions has been magnetic 
confinement of clouds (Rees 1987; Bottorff & Ferland 2000). 

However a key point is that a typical cloud needs only to ex- 
ist long enough to reprocess radiation (Rees 1987; Celotti et al. 
1992; Kuncic et al. 1996, 1997). If clouds are destroyed there- 
after and new clouds are formed, then the system can exist in 
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dynamical equilibrium even though each cloud loses its identity 
swiftly. This may apply to the broad line region. Short cloud 
survival and regeneration times can be associated with turbu- 
lence (Bottorff & Ferland 2001; Bottorff et al. 2000). The bal- 
ance between formation and destruction maintains a constant 
BLR cloud number density. 

Our analysis suggests that a nonlinear wind instability result- 
ing from the combination of launching profile + radiation pres- 
sure can lead to the formation of a turbulent biconical zone that 
might manifest itself observationally as the BLR. This biconi- 
cal structure is the region where flow streamlines cross. Such 
a bicone is remarkably similar to that inferred empirically by 
Elvis (2000) and from numerical simulations by Proga et al. 
(2000). Growing observational data corroborates the existance 
of such a structure (e.g. NGC 1068 (Arribas et al. 1996), Mrk 
3 (Ruiz etal. 2001), etc.) 

We formulate the problem in section 2.1 and describe the 
flow field in section 2.2. We analyze the linear and nonlinear 
stability of the flow and a possible scenario for thermal insta- 
bility resulting in formation of the BLR clouds in section 2.3. 
In section 3 we discuss the results and their implications for the 
AGN structure and dynamics. 

2. DISK WIND MODEL 

2.1. Formulation of the problem 

All of our calculations are performed in cylindrically sym- 
metric geometry, however we do not include axial rotation. At 
the origin we assume a black hole of mass Mbh- Fig. 1 shows 
the setup. 

We take a Keplerian, geometrically thin, optically thick ac- 
cretion disk to be located along the abscissae with the distance 
A measured along the disk. The central engine is taken to be 
a point-like source at the origin coincident with the black hole. 
We exclude disk surface illumination by the central source in- 
cluding all disk properties into the disk-wind launch velocity 
profile. 
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Our wind originates at the disk surface with launch velocity 
profile 

v z o(A) = v*(A^j , (1) 

where v* is the wind velocity at the innermost launch point of 
the disk with coordinate A* and n is an arbitrary non-positive 
number 4 . We assume that the wind propagates into a low den- 
sity medium. 

After the wind material leaves the disk, it is subjected to grav- 
ity from the black hole and radiation from the central engine. 
Proga et al. (1999, 2000) discuss one possible approach to line 
driven winds in the context of numerical modeling. To capture 
the most important qualitative features of the flow and to sim- 
plify the analysis we assume a continuum driven wind and use 
only Thomson opacities. This assumption somewhat underes- 
timates the driving force since it excludes the contribution of 
line-driving. However, on the other hand, it is an overestimate 
since typically a significant fraction of the central source spec- 
trum lies in the Compton scattering regime. To compensate for 
that we introduce a factor f T < 1, intended to decrease the total 
luminosity L. The choice of a specific value of f T depends on a 
particular situation being considered. 

In our model, a fluid parcel of density p leaves the disk with 
the velocity of equation (1), determined by the coordinate of 
its launch point, and finds itself in a centrally symmetric, con- 
servative potential field. For Thomson scattering the radiative 
force per unit volume is 

F rad (r)= 1 -^ ( auFljdv= n -?l± F ? ( 2 ) 
c H J c \r\ 

and the effective potential field is 



U(r) = U rad + U, 



grav ■ 



pGM BH I f T L 
r \ L Edd 



1 



a 

■P-, (3) 
r 



where L Edd = 1.25 • 10 38 erg/s{M BH /Mo) is the Eddington lu- 
minosity. Note, that here we give an expression for the potential 
energy per unit volume. 

2.2. Description of the flow field 

Finding a trajectory and velocity of a fluid parcel in the out- 
flow reduces to the problem of the motion of a particle in a 
centrally symmetric repulsive potential, described by (3). The 
trajectory is a hyperbola (Landau & Lifshitz 1976), and in polar 
coordinates is 

r(4>)= ^— T- ( 4 ) 



ecost 
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Here <j) g [0,7r/2]. Eccentricity of the orbit e and parameter of 
the orbit p are defined as 

2n+l 



e=l + ^, p = 2K[^ 



(5) 



where ro = Ao is the coordinate of the fluid element launch point. 
The dimensionless quantity 

K= A*(v*) 2 /2a, (6) 

is the ratio of the initial kinetic energy of a fluid particle 
launched at A* to its potential energy and measures the relative 
importance of the disk wind vs. central radiation source. (K and 
n are the key parameters in our analysis). Equation (4) can be 
considered as an implicit solution for Ao as a function of z and 
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A. For any point (A,z) we can uniquely determine the launch 
point of the streamline, Ao. Fig. 1 shows a sample streamline 
pattern. 

The velocity distribution along a streamline is 

sinh£ 



v A (A ,A) = 
v z (A ,A) = 



Ao ecosh£+l 

a , , ~~~ ~ cosh£ 

— (e 2 -l)(e+l) 

Ao ecosh£+l 



(7) 



Parameter £ e [0,+oo) is defined by cosh£ = (A/Ao)(e+ l)-e. 

2.3. Stabilty analysis 

We base our linear stability analysis for the inviscid case on 
the Fj0rtoft theorem, which is derived from the inviscid limit 
of the Orr-Sommerfeld equation, i.e. Rayleigh equation (Yih 
1969; Maslowe 1985). We have studied the outflows with the 
values of -2.5 < n < 0.0 and 0.05 <K< 2.5. In all the regimes 
studied, the disk wind is linearly stable to infinitesimal velocity 
perturbations. 

However a study of equation (4) shows that a non-linear in- 
stability can be present in a subregion of the flow. Depend- 
ing on n, a streamline inclination angle at any given point 
(r, A) is either a monotonically increasing or monotonically de- 
creasing function of Ao- The condition for the change in the 
character of monotonicity can be obtained: using equation (4) 
we can determine the inclination angle of a given streamline 
tan6(Ao,A) = dz/d\. Using 



— ( lim tan9(Ao,A)) =0, 

A— -co / 



(8) 



(/Ao V A^oo 

we find that the monotonicity changes for n = -1/2, i.e. the 
Keplerian disk is the borderline case. This means that for 
n > -1/2, each successive streamline is less inclined towards 
the abscissae axis than the preceding one and eventually inter- 
sects all of the preceding streamlines. For n < -1/2 the incli- 
nation angle is a monotonically non-increasing function of Ao 
and then the streamlines never intersect. 

The intersection of two streamlines can be approximated as 
an intersection of two jets of vanishing thickness whose inter- 
action leads to a resultant turbulent jet with an inclination angle 
lying between that of the two original jets (e.g., see (Landau & 
Lifshitz 1959)). Therefore, streamline intersection causes the 
biconical zone of compressible turbulence in the overall disk- 
wind outflow. The extent of such a zone can be estimated: since 
the streamline inclination angle is a monotonically increasing 
function of Ao, the streamline originating at the innermost ra- 
dius of the disk is the most inclined one. Therefore, it delin- 
eates the lower boundary of the turbulent biconical zone. Its 
inclination angle is 

tan6 m ,„= lim ^(A*, A) = 2VK+K 2 . (9) 

A^oo d\ 

An exact shape of the upper boundary of the bicone can be 
found by solving dz/dXo = for Ao as a function of A and 
then substituting it into the equation (4) to find the exact depen- 
dence z( A) for the upper boundary. Instead, we give an approxi- 
mate expression for the inclination angle of the outer boundary, 
found by analogy to equation (9): 



tan 0„ 



= lim A 



2(2«+l) 



(10) 
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This defines the basic geometry of the bicone: the inclina- 
tion angle 9 = (© mal +6 m! „)/2 and the divergence angle A 6 = 
Qmax-Qmin- See Figure 1. 

Typically, the bicone has a divergence angle A0 that slowly 
increases with radius. Therefore equation (10), which gives the 
limiting inclination angle of the streamline originating from the 
point Ao, depends on the choice of such a particular streamline. 
However, the increase in the divergence angle is typically small. 
Moreover, a particular streamline used to determine Q max can 
be determined independently from the properties of the accre- 
tion disk or the estimated extent of the outflow. 

Three other effects will enter in a more detailed hydrody- 
namic treatment that influence the divergence angle of the bi- 
cone. An effect which tends to make the divergence angle 
smaller than that determined using equations (9) and (10) arises 
because the resultant of two intersecting jets lies at an interme- 
diate angle between the original two. Each successive inter- 
action of resultant jets would then tend to narrow the overall 
conical structure. In addition, the change in the optical depth 
of the outflow after the formation of the turbulent bicone would 
shade the outer parts of the outflow. This will make the bicone 
narrower by increasing the inclination angle of the outer stream- 
lines due to the reduced radiation pressure. In competition with 
both of these above effects is the turbulence itself, which could 
broaden the bicone. 

The initial density fluctuations in the turbulent biconical zone 
cannot themselves represent the BELR clouds; the density and 
temperature contrasts are too small compared to the ambient 
warm highly ionized gas (WHIM in the nomenclature of Elvis 
(2000)). However, thermal instability in such a medium can 
transform initial turbulent overdensities into the BELR clouds 
(see (Burkert & Lin 2000; Hennebelle & Perault 1999) for the 
discussion of two-phase medium formation via linear thermal 
instability in optically thin regions). The formation of the com- 
pressible turbulent bicone from intersecting flows can produce 
inhomogeneities with density contrast up to 4 in the adiabatic 
case and higher if one allows for the possibility of radiative 
cooling. Such high density contrasts (compared to the linear 
density fluctuations discussed by Burkert & Lin (2000)) can 
serve as a seed for thermal instability. The instabilty proceeds 
in the nonlinear regime from the onset. Rapid cooling would 
then produce clumps. The maximum clump density reached at 
a given scale is determined by how much cooling can occur dur- 
ing an eddy turnover time on that scale, i.e. during an average 
clump survival time. 

This mechanism can lead to a dynamical equilibrium de- 
scribed by the balance between in situ formation, destruction, 
and re-formation of clumps comoving with the ambient flow. 
As long as any given clump survives long enough to reprocess 
radiation, problems of clump confinement against thermal con- 
duction and destruction via high velocity differences between 
the clumps and the ambient medium are eliminated. The dy- 
namical equilibrium of the turbulent flow, ensures a constant 
number density of BELR clouds in the turbulent biconical zone 
over observation durations. 

Finally, there is a site where a linear instability of the Kelvin- 
Helmholtz type might operate, namely the shear layer between 
the disk wind outflow and the infalling (or stationary) circum- 
polar material, that fills the cone of the outflow. Simulations 
of Proga et al. (2000) show the formation of clumpy structure 
in that boundary shear layer. They attribute it to the Kelvin- 



Helmholtz instability. However, in regimes for which the out- 
flow develops the turbulent bicone, such linear instability will 
be completely overwhelmed by the nonlinear instability in- 
duced by the intersection of streamlines. In cases when the 
bicone is not formed, the outflow density decreases from the 
maximum values near the plane of the disk to nearly the val- 
ues of the circumpolar infalling material in the upper parts of 
the outflow with the density transition occurring over the large 
range of angles. BLR clouds are not expected in this case. 5 . 

3. DISCUSSION 

We considered the linear and nonlinear stability of a disk out- 
flow and found that it is linearly stable to infinitesimal velocity 
perturbations. However, when n > -1/2 and a superedding- 
ton radiation source is present, nonlinear instability leads to the 
formation of a biconical zone of compressible turbulence. High 
density contrasts in this zone may trigger thermal instability and 
lead to the further condensation of clumps into BELR clouds. 
Turbulence plays a key role in establishing the dynamical equi- 
librium which maintains a steady cloud number density, even 
though each cloud is short lived. 

In the regimes when the outflow can develop the turbulent 
bicone, the geometry is remarkably similar to that described by 
Elvis (2000). For a disk outflow with K s=a 0.05 and an initial 
launch velocity profile slightly flatter than Keplerian, namely 
n = -0.45, the turbulent bicone inclination angle is 6~ 28° 
with the divergence angle A0 w 7°. This is in quantitative and 
qualitative agreement with Elvis (2000): (1) no absorbers are 
along the line of sight passing above the outflow; (2) along the 
line of sight passing through the turbulent bicone broad absorp- 
tion and emission line features will be observed; (3) narrow 
absorption lines (NAL) will be seen along sight lines that fall 
inside of the bicone. Figure 1 shows the geomtry of the outflow 
and the emerging bicone with the angle values found here. 

The outflow geometry in our picture depends only on dimen- 
sionless quantities K and n. However, in addition to 9 and A0 
determined above, we can also infer the inner scale by matching 
observed cloud velocities. We assume the same value of n and 
K as in the previous paragraph. We take M BH = 10 9 M Q and 
broad absorption line velocity v BAL = 10 4 km s _1 (the velocity 
observed when the line of sight passes directly along the bicone 
towards the central source). In our model vbal is the terminal 
velocity at infinity of material launched with the largest initial 
velocity and potential, i.e. at the point Ao = A*. Setting Ao = A* 
in (7) and using the expressions for va(A*,A) and v z (A*,A) in 
the limit A — > oo, we find 

A' = (11) 

V BAL 

Assuming L = 1 .5L£dd and fj = 0.7 the corresponding value of 
A* = Ao w 1000 a.u., which matches Elvis (2000). The corre- 
sponding maximum launch velocity is w 2000 km s" 1 . There- 
fore, the observed velocity in the NAL region is <~ 1000 km s" 1 . 

In principle, the paradigm presented here applies not only to 
luminous AGN but to any source with a centrally symmetric 
potential, a sufficiently luminous central radiation source, and 
a disk wind, e.g. disk winds in young stellar objects. 

We have not discussed the role of magnetic fields, disk ro- 
tation, metallicity, wind density fall off, or the physics of line 
driving. These should be considered in future work. 
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FIG. 1 . — Sample streamline pattern for a disk wind outflow. Shown is the case n = -0.45 and K = 0.06. Shaded region is the turbulent bicone formed by the 
intersection of streamlines. Values of the angles shown are discussed in the text. 



